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By solving the self-consistent system of Maxwell and density matrix equations to the first order 
with respect to nonadiabaticity, we obtain an analytical solution for the probe pulse propagation. 
The conditions for efficient storage of light are analyzed. The necessary conditions for optical 
propagation distance has been obtained. 
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I. INTRODUCTION 

Recent advances in quantum information science have 
shown that, on one hand, photons are ideal carriers of 
quantum information, and on the other hand, àtoms 
represent reliable and long-lived storage and processing 
units. In recent years quantum light storage is one of the 
cxtensively studied tasks of quantum òptics. Basic idea of 
light storage is electromagnetically induced transparency 
(EIT) Q. 

Electromagnetically induced transparency is a coher- 
ent interaction process in which a coupling làser field is 
used to made the optical dense media transparent for 
the probe field. Since its discovery, a number of new 
effects and techniques for light-matter interaction have 
appeared [2-6] . Most notably, from the point of view of 
the work presented here, particular attention has been 
devoted to ultraslow light propagation and light storage 
techniques [3-5]. 

The key concept of EIT is the dark state and popu- 
lation trapping The dark state is a specific coher- 
ent superposition state which does not contain excited 
short-living atòmic level due to destructive interference 
between two interaction paths. The dark-state is eigen- 
state of the light - atom interaction Hamiltonian, so the 
atom prepared in a dark-state can not be excited and 
cannot leave the dark-state if the interaction is adiabatic 
(Fig.l). The population trapping via applying strong 
coupling leads to the adiabatic formation of the dark- 
state. Since the interaction is realized by the light pulses 
the infuence of nonadiabatic corrections may become im- 
ortant |jj - This infuence has been studied e.g. in 
In particular, the first nonadiabatic correction con- 
nects the dark state and the bright state, so the depletion 
of the bright state, because of optical pumping, can affect 
the dark state. Despite of large amount of experimental 
and theoretical papers concerning light storage (see 
and citations there), and applications in quantum infor- 
mation science, the influence of decoherence level width 
on information carried by stored light is studied insuffi- 
ciently. 

In this work we present theoretical study which dis- 
cusses and explains influence of all relaxations on probe 
propagation both analytically and numerically (it is es- 



sential in especially, sòlid state systems [Ï0|1. The goal 
is to study comprehensively how the depletion of bright 
state will affect the pulse propagation in an EIT media 
and light storage in particular. By solving the couplcd 
system of Maxwell and density matrix equations to the 
first order of the nonstationary pertrubation theory with 
respect to nonadiabaticity and decoherence we obtain an- 
alytical solution which completely describes the probe 
pulse propagation and is consistent with the recent light 
storage experiments. 

The paper is organized as follows. In section II the 
bàsic equations are written down and the probe pulse 
propagation equation is derived and analyzed. In section 
III and Appendixes the analytical solution for counterin- 
tuitive pulse switching order and for matched pulses are 
obtained and the asymptotic solutions discussed. Sec- 
tion IV deals with the physical consequences of the ob- 
tained solution, namely the necessary conditions of the 
pulse storage and retrieving, also the numerical results 
are demonstrated. In section V we consider the trans- 
verse relaxation of the coherence induced in the medium. 
Section VI concludes the paper. 



II. BASIC EQUATIONS 

Figure 1 shows a schematic diagram of the atòmic sys- 
tem in the EIT basis: media of three level àtoms interact- 
ing with two làser pulses E p = A p cos (u> p t — k p z + ip p j 
(probe) and E c = A c cos (u> c t — k c z + tp c ) (coupling). 
The probe field resonantly conneets the state |1) to the 
state |3) and the coupling field conneets |2) to |3). The 
Hamiltonian of the system in the rotating wave approxi- 
mation is: 

H = hA<r 33 — hfl p <j 31 — hQ p a 32 + H.c, 



where fL 



A, 



"p.c — z — (* — 1, 2) are the respective Rabi 

freqüències, er^ = are the atòmic transition oper- 

ators, A = lü p — uj 3 i — lo c — lo 3 2 is the detuning of the 
pulse freqüències from the upper level and (i = 1, 2) 
are the dipole moments of corresponding transitions. 

We assume that: (i) the probe field is weak as com- 
parec! to the coupling pulse field íl p << O c ; (ii) the in- 
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teraction is adiabatic (O c T >> 1, where T is the interac- 
tion duration). Then, the atòmic density matrix equation 
may be written as 



Pzi = -rp 31 +íO p + ffl c /3 21 , 

p21 = ïíïcPSD 

P11 = !» 

P22 = P33 = P32 = 0; 



(1) 



where T is the width of the upper level which is the sum of 
the spontaneous decay and transverse relaxations rates. 
It is supposed that interaction is fast enough to neglect 
the decoherence between metastable levels (sec. III, IV), 
or to take it into account to the first order (sec. V). 

The propagation of the pulses is governed by the 
Maxwell equation for slowly varying amplitudes, 



where qi 
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, N is the atòmic number density. 



System of equations (yi can be reduced to one equation 
for p 31 



P31 



■Paí 



P31 



Paiçf) = iÇl c' e ( 3 ) 



where 9 — — — is the commou used notation for the so 

called mixing angle. Influcncc of the first two terms in 
(PJ) can be neglected if we confino to only first terms with 
respect to the nonadiabaticity (i.e. (fi c T)~ << 1 is 
neglected) . Influence of the fourth term in © is essential 
parameter only under the assumption 



TT » 1. 



(4) 



The meaning of the condition is obvious: under the 
condition of complete adiabaticity relaxation does not af- 
fect the pulse propagation (dark-state), but taking into 
account first nonadiabatic correction, has essential influ- 
ence. The relaxation can be neglected when TT < 1. 
Finally, by substituting the Maxwell equation @ into 
one gets pulse propagation equation in wave variables 

X — X ^ í - — í x I ' C-) 



1 d 2 9 
T\dxdt 



q L· de ae = 
dt dx 



fi 



(5) 



where notation Ti = — - is used. In this conncc- 

. r. 

tion the coherence dynamics is governed by the following 
equation: 



P21 = - r iP2i - 



(6) 



Thus Ti is the coherence decay rate, or width of EIT 
resonance, due to applied coupling. Under the condition 



T\T » 1 



(7) 



equation © has the well known quasi-stationary solution 
P21 = —9 p|, and the equation © passes to the dark 
state polariton propagation equation. The condition (Q 
means, that width of EIT resonance exceeds the spectral 
width of the probe. 



III. SOLUTION OF PROPAGATION EQUATION 

The obtained probe pulse propagation equation |5j is 
solved by the method presented in Since (JHJ) is lin- 
ear in 9 and Q c (t) is independent of x, it can be solved 
by using the Laplace transform with respect to x. The 
solution of JSJ for 0's Laplace image can be found easily: 



(M) 



9p + ri6>ç 
s + q p /T 



B (s, t, íx) dti + c (s) B (s, í, —00) 

(8) 

/ t \ 

J Tidt' , c(s) is an 



where B (s. t, U ) = exp , 

integration constant that is determined by the initial con- 
dition, c (s) = 9 (s, —00). If pulses are switched in coun- 
tcrintuitivc sequence (coupling turns on earlier than the 
probe does) then c (s) = 0, since 9 (z, —00) = 9q (—00) = 
(see appendix A). 

Space time evolution of the probe pulse is obtained by 
implementing the reverse Laplace transform in J5J ■ 



>{z,t) 



(0o (í) r x + Ò (í)) x 



(9) 



x exp (-z - a (h,t)) I (2yJ za (h,t)j , 



where z = is propagation distancc normalized to 

linear absorption factor and, for convenience, the nota- 

í 

tion a(ti,t) = j Ti (t')dt' is used. By using the con- 
ti 

dition za(ti,t) » 1 one can substitute the modified 
Bessel function by its asymptote, so the solution JÜJ re- 
duces to the following: 



'OM) 



x exp 



dh (flo(í)ri (*i) + 0o(í) x 



(10) 



(- (y/z- V / o(ti7í)) 2 J («a (íi,*)) 1/4 . 



Depending on optical propagation distance z two sim- 
ple asymptotes for (|10f> can be obtained (see appendix 
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B). The first is the case where 



(11) 



Ti m is the maximal value of Ti (t) (see also 0)- Un- 
der condition , solution © reduces to the dark-state 
polariton propagation solution with correction in Q: 



e(z,t) = e (0 



i 



(0 



(12) 



where Ç is the non-linear time determined by 
t 

J Vl 2 c (ti) dti — q p x |íj. Note, that turning off the cou- 

í 

pling íl c (t) does not reduce Ti (£) to zero, since £ retards 
from t. In this case, as it will be shown below, the infor- 
mation stored in the medium can be well retrieved (see 
sec. IV). 

In the case of condition reversed to (|llfl . 



the solution reduces to 



(z,t) = iïexp (- (yz- ^qí(ío,í)^ j(za(t ,t)) 



(13) 



,-1/4 



(14) 

where to is the maximal value of 9$ (t) , R = 

oo 

J Ti (t') 9o (f) dt' and does not depend on time. We 

— oo 

emphasize that for propagation distances meeting the 
condition i|13|) . the obtained pulse loses all the informa- 
tion about its initial temporal shape, since the right hand 
side in l(T4^l does not contain time dependent 9q- 

IV. DISCUSSION 

In this section the probe pulse propagation dynamics 
obtained from the analytical solution (JÜJ) is presented. 
First of all we consider the case of constant coupling field. 
Shape of the initial pulse is chosen to be double-humped 
in order to visualize the propagation dynamics. For the 
small propagation distances when the condition (|ll|l is 
met influence of V is negligible (Fig.2a). When the con- 
dition (fTT]) is violated, the influence of upper level width 
becomes essential as one can see from Figs. 2b, c. Thus 
influence of T breaks the adiabaton propagation regime. 

As it was mentioned above, propagation over very long 
distances (|13[) leads to the lost of the information on the 
initial pulse temporal shape. This can be seen in Fig 3, 
where propagation over the same distance of two pulses 
with different temporal shapes but with the same initial 
area is depicted. By propagating over very long distance 
(|13|l they lose any information about their initial tempo- 
ral shapes. 



In Fig. 4 we show that the increase of Ti suppresses 
the smearing of the probe. This is caused by the decrease 
of the bright state population and hence leads to the de- 
crease the influence of T on pulse propagation. Note, that 
in the dark-state propagation regime the pulse temporal 
shape does not depend on coupling field amplitude or on 
unstable level width. 

Summarizing presented results one can see that to min- 
imize the pulse smearing during its propagation one has 
to either increase Ti or decrease the propagation distance 
z. The situation changes dramatically for the light stor- 
age and retrieving process (fl c ^ const). 

It is known, that pulse can be completely stored and 
retrieved from the medium if the medium length and Ti 
meet the condition (see for example 0): 



(15) 



For efncient storage and retrieving the condition 111|) 
also has to be met. Combining this two nonequalities one 
gets that to completely store and well retrieve the light 
pulse, T\ has to meet the following condition: 



T lm T » 161n2 >> 1. 



(16) 



Therefore, influence of the second term in H12f) is in- 
sufficient when the condition IjltjD is met. 

Storage and retrieving of thc light pulse for different 
propagation distances under the condition Ijltjll is de- 
picted in Fig 5. For small propagation distances when the 
condition (|15|l is violated only the falling edge of the pulse 
is stored and can be retrieved (Fig. 5a), because when 
this edge enters the medium, the leading edge emerges 
already. For larger propagation distances when the con- 
dition (|T5)l is satisfied the whole pulse can be stored and 
then well retrieved by turning on the coupling field. 

Let us now consider the case when the condition (|lfc>|> 
is not met (Fig. 6). For small propagation distances 
when the condition l|ll|) is satisfied but l|15|) is not, only 
the falling edge of the pulse can be stored and retrieved. 
Propagation over longer distances brings to satisfying of 
(|15fl and violation of (|llfl . Thus the whole pulse can be 
stored but the retrieved pulse temporal shape is smeared. 

We present finally comparison of the experimental re- 
sults with our analytical solution. In Fig 7a the experi- 
mental data of storage and retrieving of light pulse from 
[g are presented. Curve in Fig 7b is plotted from our 
analytical solution JjJ}: all parameters correspond to the 
conditions of the experiment. One can see good con- 
sistency between experimental data and our analytical 
solution (Note that the storage in case of experiment () 
is incomplete as was discussed above). 



CONSIDERATION OF p 21 TRANSVERSE 
DECAY 



In this sections we take into acount the quantity 7T in 
first order. This leads, instead of (Q, to the equations. 
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p 31 = -Tp 31 +iCl p + iiï c p 21 , 
Px\ = !» 

Pl2 = P33 = P32 = 0; 

Thus, probe pulse propagation equation is written as 
follows: 



1 d 2 9 89 q p 09 q pl 



Tx dxdt dx T x (r + 7) dt T í {V + 7) 
where Ti (t) now is 



0, (17) 



+ 7 r + 



Ti (t) 



r + 7 



(18) 



As results from (|18|l , to completely stop the light in the 
medium (l'i = 0) one should turn off the coupling field 

ü c at the rate T (i.e., r + ^ = 0). 

By performing the stated above Laplace transform pro- 
cedure one obtains analytical solution of the equation 
(|17|l in the form 



6{z,t) = J dh [Oo + rtfo) x (19) 

— OO 

x exp I -z - J Txdt' J J (2y/za (íi,í)) 



where a (*i,í) = / Ti (f) - jdt' and z = 

More detailed analysis of the expression l|17(l will be 
performed in a subsequent publication. 



VI. CONCLUSION 

We have considered the propagation, storage and re- 
trieving of the light pulse in EIT media by taking into 
account all dephasing rates. From coupled system of 
Maxwell and density matrix equations we derive the 
probe pulse propagation equation, which in particular 
case passes into the dark-state polariton propagation 
equation. We find an analytical solution and analyzed 
its physical consequences. We derived a simple asymp- 
totes of the solution, and showed strong dependence of 
light pulse temporal shape on optical propagation dis- 
tance in the presence of relaxations. We demonstrated 
that an efficient storage of light is possible by choosing 



appropriate coupling intensities and optical propagation 
distances. Finally, we compared our solution with exper- 
imental data and showed that our solution is well consis- 
tent with the recent experiments. 
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APPENDIX A: 

To find the c (s) in general case we examine the equa- 
tion ©. 

. .30 iq„ , . 

From this equation one obtains — = — p 3í (x, t), 

ox ll c 

and the solution can formally be written as 9 (x, —00) = 
iq 

ü 



J p 31 (x', —00) dx' + Qq (—00) . This can be simplified 

c 

by taking into account that there is no dipole moment in- 
duced before the interaction is turned on (p 31 (x, —00) = 
0), 

9 (x, —00) = 9q (—00) 

For the counterintuitive order of pulse switching we 
have 6*o (—00) = 0, hence c(s) = 0. In the case of si- 
multaneous pulse switching 9q (—00) = const, therefore 

. . const 

c(s) = . 

s 

In the case 0q = const from general solution © one 
easily obtains the so-called matched pulse propagation 
regime 9 = 8 [Ü- 



APPENDIX B: 

The integrand function in l|10|) is the product of two 
time dependent bell-shaped functions. This kind of inte- 
gral is evaluated easily when the temporal width of one 
function is much larger than the width of the other one. 
Then by the first order saddle-point technique the narrow 
function can be approximated by delta function. 

The time widths of the first and second functions in 
(|10fl are, respectively, the interaction duration T and 
4V7In72) 

Tlm 

On the one hand, if the condition l|ll|l is met, the inte- 
gration of (|lü|) gives the well known dark-state polariton 
solution l|12|l . On the other hand, if the condition l|13l) 
is satisfied, integration of i|10fl by the saddle-point tech- 
nique brings to the asymptote (|14|l . 
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Fig.l Schematic diagram of the three-level atòmic sys- 
tem in the basis of EIT in the dressed-state (a) and the 
bare-state (b) representations. The field íl c couples the 
\bright) state to the upper level |3), and another field 

with Rabi frequency 9 couples the \dark) and \bright) 
states. 

Fig.2 Probe pulse temporal shape after propa- 
gation over different propagation distances in the 



case of constant coupling field; T\T — 40. Ini- 
tial temporal shape of the probe has been cho- 
sen as Ü p /Y = 0.012 exp (- (Tí/100 - 7.5) 2 ) + 

0.01 exp (- (rt/100 - 10) 2 ) . For small propagation dis- 
tances when the condition is met the influence of 
upper level width is small (a). Violation of the condition 
(fTT|) increases the role of the upper level width: (b), (c). 
Here and below the dotted curve is the initial temporal 
shape of the probe. 

Fig.3 Propagation of two probe pulses with differ- 
ent temporal shapes but with the same area over very 
long distances l|13|l in the case of constant coupling 
T\T = 4. One can see that during propagation pulse 
loses all information about its initial temporal shape. 
The shape of the probe in the case (a) is the same as 
in Fig.2, in the case (b) it has been chosen as Íl p /T = 

0.022 exp (- (Ti/lOO - 8.5) 2 ) . 

Fig.4 Probe pulse temporal shape after propagation 
over the same optical distance z = 20 for different con- 
stant coupling amplitudes. Increasing of coupling ampli- 
tude decreases the influence of upper level width. The 
shape of the probe is the same as in Fig.2. 

Fig.5 Storage and retrieving of the probe for differ- 
ent propagation distances, when the condition (|16|) is 
met. The coupling field is chosen as íi c /r = / (i) 
(a), where / (i) = 1 for Tt < 1000 and Tt > 2500; 

/ (í) = exp (- (Tí/100 - 10) 2 ) +exp (- (Tí/lOO - 20) 2 ) 
for 1000 ^ Tt ^ 2500. For small optical propagation dis- 
tances only the falling edge of the pulse is stored and 
retrieved (b). For larger z the whole pulse can be stored 
and well retrieved (c). The shape of the probe is the 
same as in Fig.2. 

Fig.6 Storage and retrieving of the probe for differ- 
ent propagation distances when the condition i|16|) is vi- 
olated. The coupling field has been chosen as íí c /r = 
/ (í) / VTÜ, T lm T = 40 (a). Either only the falling edge 
can be stored and well retrieved (b) or the whole pulse 
is stored but loses its initial shape (c) . The shape of the 
probe is the same as in Fig.2. 

Fig.7 Experimental data from [6|, (a) and our analyt- 
ical solution for parameters corresponding to the exper- 
imental conditions (b). ft p /T = exp í- (rí/400 - 10) 2 ) 

for Tt < 1000, and Q p /T = exp (- (3rt/100- 10) 2 ) for 

rt > íooo. n c /r = 0.1/ (rt) , r lm T = 5, time axes is 

plotted in [i sec . 
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